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Abstract 

Rapidly convergent expansions of a one-loop contribution to the partition function of quantum fields with 
ellipsoid constant-energy surface dispersion law are derived. The omega-potential is naturally decomposed 
into three parts: the quasiclassical contribution, the contribution from the branch cut of the dispersion law, 
and the oscillating part. The low- and high-temperature expansions of the quasiclassical part are obtained. 
An explicit expression and a relation of the contribution from the cut with the Casimir term and vacuum 
energy are established. The oscillating part is represented in the form of the Chowla-Selberg expansion for 
the Epstein zeta function. Various resummations of this expansion are considered. The developed general 
procedure is applied to two models: massless particles in a box both at zero and non-zero chemical potential; 
electrons in a thin metal film. The rapidly convergent expansions of the partition function and average 
particle number are obtained for these models. In particular, the oscillations of the chemical potential of 
conduction electrons in graphene and a thin metal film due to a variation of sizes of the crystal are described. 
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I. INTRODUCTION 



Evaluation of one-loop corrections is a cornerstone mean allowing to describe the quantum 
corrections to physical properties of many-particle systems. In many cases, the one-loop corrections 
together with their renormalization group improvements provide the only attainable information 
about quantum features of the system which can be analyzed analytically. An identity of the formal 
expressions for these corrections makes to use the same computational techniques in different fields 
of physics from astrophysics and cosmology with characteristic energy scales up to hundreds of 
GeVs down to condensed matter physics with energies of the order of a few eVs and lesser. The 
calculations of the one-loop corrections are confronted with certain mathematical difficulties which 
are common for all these branches of physics. The problem becomes extremely complicated when we 
investigate quantum fields with imposed non-trivial boundary conditions and/or interacting with 
background fields of a complex structure. 

There is a well elaborated quasiclassical approach known as the heat kernel expansion for how 
to calculate the one-loop corrections both at zero (see for review ||l|) and finite temperatures and 
chemical potentials |||, for almost any background fields and boundary conditions. Nevertheless, 
it possesses some weaknesses. This method is mainly developed for particles with a relativistic 
dispersion law (see, however, |^). Besides, and what even is more important, it does not catch 
exponentially suppressed terms which are substantial and observable at certain values of external 
fields. In condense matter physics, the classical example of these contributions provide the oscillat- 
ing terms in the ri-potential of conduction electrons in metals. It was apparently Landau |6| who 
first observed that an applied magnetic field may result in the oscillations of the partition function 
of conduction electrons. These oscillation were experimentally discovered by Schubnikov and de 
Haas in conductivity and by de Haas and van Alphen |^ in magnetic susceptibility. Later on, a 
quasiclassical theory of these oscillations was developed for the Fermi liquid quasiparticles with an 
arbitrary dispersion law p|~]T4||. 

In this paper, we are going to consider the free quantum fields subjected to certain boundary 
conditions. The interaction with an environment and the self-interaction are taken into account 
only in the form of the dispersion law which we choose to be 

E{p) = Lo{g'^{p^ + ai){pj + aj) + m^), (1) 

where gij^ Ui, and m are some constant quantities, gij being symmetric and positive definite, and 
uj{x) is some smooth function increasing at a large argument. Henceforth, the Einstein summation 
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rule is assumed and indices are risen and lowered by the "metric" g'^^ and its inverse. This metric 
is related to the effective mass tensor, 

- -'W ^ ^ 

when the latter is defined. Such a problem definition is typical for the Fermi liquid theory of 
conduction electrons in metals where the dispersion laws take exotic forms largely deviating from the 
familiar non-relativistic and relativistic ones (see, e.g., |[Tl| ). Assuming relation (|l]), we consider the 
simplest modification of the dispersion law such that its constant-energy surface, and, in particular, 
the Fermi surface, has an ellipsoid form. Of course, the form of the dispersion law (Q) chosen by 
us is not accidental. Only in this case, we can reduce the problem of analytic evaluation of the 
partition function to a one-dimensional one. This, in turn, admits us to investigate characteristic 
features of the fi-potential in details analytically going beyond the quasiclassical approximation of 
the general theory of oscillations ||T^ . Our aim is to derive the rapidly convergent expansions with 
a controllable error of the one-loop partition function for particles with the dispersion law 
In calculating the il-potential, we shall see that it naturally falls into three pieces 

r2 = -)- rig, ^iq = -i- ^os- (3) 

Here is the quasiclassical (main) contribution. A discreteness of quantum numbers is completely 
neglected in this term. It is that contribution which is described by the naive heat-kernel expansion. 
The term Qq is an essentially quantum contribution. It is decomposed on the contribution from the 
branch cut of the dispersion law (|l]) (if exists) and the oscillating term which is given by the sum 
over the Matsubara frequencies. The latter term oscillates with the chemical potential, while the 
former does not. Notice that, in this paper, we do not regard the problem of a vacuum energy and 
the general expression (^) for the J7-potential is assumed not to include the divergent vacuum term. 
Nevertheless, as we shall see, the vacuum contribution can be found from the high-temperature 
expansion of (^), and the Casimir term |15| is contained in the contribution 



In Sec. H, we develop a general procedure for obtaining of rapidly convergent expansions of 
the partition function. We start with the definition of the fi-potential for particles (boson or 
fermions) with the dispersion law (|l|) and, by use of the Poisson summation formula, bring it to 
the form (^). Then, every term in the representation (|3|) is studied separately. The oscillating 



term J7os reduces to a sum of the Epstein ^-functions [ |lq | and, in Sec. [I A , we, in fact, rederive 



the Chowla-Selberg expansion |17-26| for it. Contrary to the standard ("-function approach to 
evaluation of the one-loop corrections, we do not assume the CPT-symmetry (an identity of the 



particle and anti-particle dispersion laws) in our calculations. In this section, we also single out the 
contribution from the cut and establish its relation to the Casimir term and the vacuum energy. In 
Sec. p A 1| , we represent the quasiclassical contribution 0,q in the form analogous to that for the 
essentially quantum contribution 0,q, but assuming an analytical continuation in the parameter d, 
which initially characterizes the space dimension. In the next Sec. [I B| , we analyze an effectiveness 



of the obtained expansion and derive the restrictions on the spectrum parameters under which 
the expansion converges rapidly. These conditions also allow us to find a region in the parameter 
space where the oscillations become significant. In order to extend a range of applicability of the 



expansion, we further study its various resummations. Sec. II C is devoted to this problem. Here we 
consider three cases: i) the poles corresponding to the Matsubara frequencies merge and constitute 
a cut; ii) the problem reduces effectively to a one-dimensional one like, for example, in the case of a 
thin film with macroscopic transversal dimensions; iii) one or several poles approach the real axis so 
that their contributions to the partition function are not exponentially suppressed. In Sec. II D, we 
turn back to the quasiclassical contribution and derive its low- and high-temperature expansions. 
In particular, the obtained high-temperature expansion generalizes the known one ^ for the 
relativistic dispersion law to the dispersion law of the form (^. Here we also explicitly write out 
the finite and logarithmically divergent contributions to the one-loop vacuum energy for particles 
with the dispersion law (|^). It turns out that the logarithmic divergencies which are important for 
the renormalization group method arise only for certain dispersion laws, the relativistic spectrum 
being one of them. 

In Sec. m, we apply the general formalism to two models: massless relativistic particles (a 
linear dispersion law) in a box Sec. pIAl , and non-relativistic particles in a thin metal film Sec. 



IIIB. The latter model is a classical subject of the Fermi liquid theory of conduction electrons 



in metals and was studied in many papers and books (some of them are |, IH-Hl). As it 

was mentioned, our purpose in this case is to obtain a rapidly convergent expansion more accurate 
than the quasiclassical result of the general theory of oscillations Here we obtain the closed 

simple expressions for the il-potential and the average number of conduction electrons which are 



valid with an exponential accuracy. Also in Sec. HI B , we consider a dependence of the chemical 
potential on the thickness of the metal film. It is known 11, 12| that the partition function 
oscillates when one changes the sizes of the system. We derive the simple expressions for a period 
and amplitude of these oscillations. The first model is related, for example, to photons confined 



to the ideal metal rectangular parallelepiped |15, when the chemical potential vanishes, and 



to electrons near the Dirac points ||3^, ^] in graphene (graphite) when the chemical potential is 



not zero. Despite a study of these subjects also has a long history, they remain hot nowadays. 
The investigation of thermodynamic properties of massless particles in pistons is provoked by the 
experimental observations of the Casimir force (see for review ||34|), while the second subject is 
concerned with graphene [^^-08|. 



In considering the conduction electrons in graphene and a thin metal film, we make certain 
approximations. First, we completely ignore an interaction of the electrons with the electromagnetic 
field not produced by the core electrons and nuclei, the latter being taken into account only in 
the form of the dispersion law. This is a standard approximation. A theory which includes this 



interaction both at zero (see, e.g., |3£, |42|) and finite temperature |41| is known. In connection with 
the Casimir effect, its path-integral formulation was recently given in [42| in the Dirac approximation 
(relativistic dispersion law). Second, describing the oscillations of the chemical potential due to a 
variation of sizes of a graphene specimen, we suppose that the additional contributions to it resulting 



from a deformation of the lattice |37, 43, 44 1 are approximately constant. This changes are rather 
small for the crystal sizes and deformations we study. Moreover, we can simply add them to 
the electronic chemical potential relaxing thereby the constancy assumption. Some experimental 
investigations of the graphene properties under pressure where the deformations of its lattice are 



realized can be found in |45|. We make the same assumption in studying the oscillations in a thin 
metal film. Here this approximation is a standard one since a variation of the dispersion law of 
conduction electrons is small under reasonable deformations of the specimen. Third, we do not 
include in our consideration any disorder effects supposing that the considered crystals are close 
to ideal. In a certain approximation, these effects are taken into account by a convolution of the 
obtained expressions for the il-potential over the chemical potential with the Lorentz distribution 



|41|. 



It will be convenient to use the following notation. We denote by /i and call it the dimensionless 
chemical potential the quantity /3/i, where /i is a customary chemical potential. Besides, the Q- 
potential is defined as 



:=lnZ(/3,^,e), Z:=Spe-^^+^^, (4) 



where H is the many-particle Hamiltonian, is a conserved charge and ^ some other parameters 
characterizing the spectrum of the Hamiltonian. 



5 



II. ONE-LOOP SI-POTENTIAL 



In this section, we derive various rapidly convergent expansions of the series 

Id = J];in(l + e^-^(P)), pGZ^ (5) 
p 

with the dispersion law (||). The sum in Eq. (^) is taken over all the integer pi, and, for convenience, 
we define the dimensionless chemical potential fj, in such a way that E{p) > for any p. In the case 
we study, the one-loop contribution to the il-potential is proportional to this series or expressed in 
terms of the series of such a type. Usually, the series (^) corresponds to the partition function of 
quantum fields subjected to the periodic boundary conditions. As for bosons, the logarithm of the 
partition function differs from Eq. (||) by the overall sign and by the sign at the exponent in the 
logarithm. The latter sign can be simply earned by the replacement — )• fi — iir. Making use of the 
representation of the 5-function in terms of the Fourier series, we can write (the Poisson formula) 

Id = Y. [ <fpe~'^^'P^'i' ln(l + e^--^(P)) = e^^*'''?' ( dV"^^*^'''' + e^-'^(p'+™')), (6) 
where g is a multiindex taking the same values as the multiindex p. 

A. Expansion of the oscillating term 

Let us consider in detail the integral 

Il = j d'^pe-'^™' ln(l + e/^-'^(p'+™')) =: j dV^^^'^'^W)- (7) 

At first, we reduce it to the one-dimensional integral 

ll = 2Ti^ dppja/2~,{q,p)fip^), (8) 
Jo 

where g := det Qij, p := \/p^ and q := \fcp' . Hereinafter, we use a convenient notation for the 
expressions containing the Bessel function Z^(x) 

Zuiq,p) := (-) Z,{2TTqp). (9) 



.1. 

In particular, these functions obey the useful recurrence relations 

= 2vr..-r(.,p), = -2vr..,,(,,p). (10) 

pop qoq 



The integral is an analytic function of the complex variable d at Red > 0. Besides, we shall 
assume that 

lim = lim d,J^, = 0. (11) 

This conditions will allow us to restore the function by its second derivative. The analytic 
continuation of the function to the whole complex d-plane can be achieved by the standard trick. 
We pass from the integral over the real semi-axis to the integral over the Hankel contour 



7 / dsj,,/2-iiq,s'/^)f{s), (12) 



where the contour H runs from +oo a little bit higher the real axis, encircles the origin and then 
goes to +00 a little bit lower the real axis. When d is an even number, we have to take the limit 
in Eq. (|l2|). We also assume that we can choose such the branch of the function f{s) which has 
not a cut on the real semi-axis. For any reasonable a;(s) the integral ( p!2| ) converges uniformly with 
respect to /x < ^uq < 0. Then we see from the representation (^) that the "boundary" conditions 



(11) are satisfied for any d G C. 
Consider the function 



f{p') ■■= dlfip') (13) 



as an analytic function of the variable s := in the complex plane. We shall suppose that the 
function Ijj{s + m?) defining the dispersion law is analytic on the whole s-plane or possesses the 
branch point Sc < on the real axis. A generalization to the case when the function uj{s + m?) 
owns more than one branch points will be obvious. It is useful to choose such a branch of a;(s + m^) 
which respects the Schwarz's symmetry principle 

uj*{s) =io{s*). (14) 

Then the function f{s) also has this property. The cut picking out such a branch of u}{s) goes along 
the real axis from — oo to Sc- In addition to the mentioned branch point, the function /(s) has 
singularities (the poles of the second order) in the points 

u}{sk + m^) = fi + iuk, cjfc = 7r(2/c + 1), A; G Z, (15) 

where u}{s) denotes the branch chosen by us and Wfc are the Matsubara frequencies. The solutions 
Sk of Eq. (ITsI ) can coincide at different k. Moreover, Eq. ([Tsl ) can have a non-unique solution at 
the fixed k or have not any one. Due to the symmetry (^4|), the poles appear by the complex 
conjugate pairs. Further, we shall denote by all the different solutions to Eq. (ffsl). 



Let C be the contour enclosing the poles (15) in a counterclockwise direction and Coo be the 
circle with a radius tending to infinity. Then, by the use of the Cauchy formula, /(p^) can be 
represented in a form analogous to the Kallen-Lehmann representation for the Fourier transform of 
an exact propagator 



c 



ds f{s) 
27rip2 — 



ds f{s) 

2711 p2 — 



ds f{s + ie) — f{s — ie) 



2m 



— s 



c 



ds f{s) 
27ri — s 



(16) 



where the contour C encloses only the singular point s = p'^, and e is an infinitesimally small 
positive quantity. Also we assume that f{s) is integrable in a vicinity of the branch point s = Sc- If 
the function uj{s) has no branch points, the second integral should be omitted. If f{s) tends to zero 
almost everywhere at \s\ — )• c« then the last integral is also absent. In particular, the last property 
is fulfilled for the function (13) if 



lim \uj{s)\ = oo. 

|s|^oo 



(17) 

Otherwise, we need to make subtractions in Eq. (Il6|). Discarding the integral over the large circle 
in (p!6l), we arrive at 



Sc 

f°° f ds f{s) f 

dlll = 27.^gj^ dppj,/,_^{q,p) 

IC' -oo 



ds Im/(s + ie) 



TT 



p"^ — S 



dsf{s) — 2i / dslm.f{s + ie) 



LC" 







jd/2-i{q,p) 



The change of an integration order in the last equality is valid when Red G (0,3) at g 7^ 0, and 
Re d G (0, 2) at (7 = 0. If the physical dimension does not fall into these intervals then we shall make 
all the calculations as it would be in these intervals and set its value to the real one in the final 
result only. These operations are legitimate by virtue of a uniqueness of an analytic continuation. 

Further, we shall consider the case (7 7^ 0. Notice that the integral factorized in (^) can 
be interpreted as the inverse Fourier transform of the free propagator with the mass s^/^ in a 
d-dimensional Euclidean space. Now we apply the formula ([EH|, p. 434) 



> 0, Rei/ e -1,2a + 



(19) 



(p2-s)"+i 2r(a + l) """^"^'^ " ' ~V ''' ' 2^ 

and s does not lie on the real positive semi-axis. The cut of the function s^/^ is also along the real 
positive semi-axis so s^/^ takes the values in the upper half-plane. At a = we obtain 



f,2rq _ Vd 
Vd- 2 



00 

j dsf{s)dsh^^]^{q,s"^)-A j ds\Tnf{-s + ie)kd/2-i{q,s^'^) 



(20) 



c 



where we have used the recurrence relation (10) and the relation of the Hankel function with the 
Macdonald one 



i7r/iW(g,ip) = 2k„{q,p). 



It also follows from Eq. ( |T9| ) that 



(21) 



(22) 



when v and q are real. 

Now we take into account an explicit form of the function f{s). Then, for any function ip{s) 
analytic in a neighbourhood of the point Sfc, the following relation holds 



res [f{s)ip{s)] 

S = Sk 



dJdaSkfiSk)], 



(23) 



where are specified by Eq. (15). Bearing in mind the boundary conditions (11), the integral /| 
is rewritten as the sum of contributions from the residues and the cut 
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OO 

y c^s^d/2-i(9,'5^/^)Im 



In cosh 



2 a;+(s) - ii 



ITT 



Y.^^dU'^^Pk] 



oo 



ds 



^-*"(s)fcd/2-i(9,s'/') + ^^tanh "^"'^".^~^ fc, 



TT 



d/2 



iq,s 



1/2^ 



i, ' 



(24) 



where we define pk := s^/"^ and the function w^(s) := uj{'m? + ie — s). The summation is taken 
over all the residues (|l5|) of the function f{s). According to our definition of the square root, 
Impfc > and, consequently, the contributions from the poles and the cut to the 0-potential are 
exponentially suppressed. In the last expression for the variable d is supposed to be real. In 
this expression, we explicitly singled out the Casimir contribution from the cut to the logarithm of 
the partition function (the first term in the square brackets) and also integrated the reminder by 
parts assuming that the discontinuity of the integrand goes to zero at the branch point. A bosonic 



variant of formula (24) is produced by the substitution ^ — t- ^ — ivr and change of the overall sign. 

Some comments on the expansion ( p4[ ) are in order. If the function cj(s) takes the values on the 
opposite banks of the cut which differ by a sign only, as, for example, for the relativistic dispersion 



law, then an imaginary part of the second term in the square brackets in the expansion (24) becomes 
an odd function of fi. Therefore, if the considered model allows for the antiparticles and we include 
their contribution to the Q-potential, these terms cancel out and the contribution from the cut is 



given solely by the Casimir term. For certain dispersion laws and values of the chemical potential, 
the situation may occur when some poles in the s-plane appear precisely on the cut of the function 
u){s + m?). In this case, the integrals in (p^) are understood in the sense of principal value, while 



the contribution from these poles enters the expansion (24) with the factor 1/2. 

Notice that, in order to obtain the contribution from the poles to the sum (^), we could use the 
well-known representation of the Epstein C-function in terms of a series in the Macdonald functions 
P^7|-p6|. With this observation, we would come to the expansion (24) without the contribution from 
the cut at one stroke. Also, if the function u}{s) has no branch points, we could use the general 



formula |47| for integrals of the type (|8|) and obtain the expansion (p^). 

Let us estimate a behaviour of the different contributions to ( p^ ) at high temperatures (small 
/?). In the fermionic case, the contributions from the poles tends exponentially to zero for the 
dispersion laws uj{s) with the power-like asymptotics at high energies. This holds for bosons as well 
excepting the contribution from the pole with zero Matsubara frequency. At the fixed chemical 
potential /i = its contribution to the average energy is proportional to /3~^. In the fermionic 
case, the contribution of the second term in the square brackets tends linearly to zero at /3 — )• 
and fixed jl. For bosons, the contribution of this term to the average energy behaves like /3~^ at 
the fixed chemical potential fl. The first term in the square brackets (the Casimir term) depends 
linearly on the temperature. So, its contribution to the average energy is independent of /3. It is the 
Casimir term which determines the leading asymptotics in /3 of the essentially quantum (i.e. non- 
quasiclassical) thermal correction to the average energy in the fermionic case. As far as the bosons 
are concerned, this leading asymptotics is determined by the above mentioned terms proportional 
to At low temperatures and E{p) > 0, all these terms are absorbed by other contributions to 
the average energy so that, in sum, it tends exponentially to zero. 

The Casimir term is exactly canceled by the analogous term coming from the vacuum fluctu- 
ations. It is easy to understand from the following general observation. An expression for the 
one-loop vacuum energy can be obtained from the high-temperature expansion for fermions at the 
vanishing chemical potential since 

Efi/3,0) = - d,mj{p,^.))\^^, = y: ^^^^ ^„ E (25) 

n n 

where En is the energy of the mode n, the index / marking the statistics. In this method for 
obtaining the vacuum energy, the Fermi-Dirac distribution plays the role of a regulator resulting in 
the energy cutoff. Let the contribution to the average energy at zeroth power of /3 takes the form 

Ef{/3, fi) = a[J(/x) + a?(/.) ln(/36?(/.)) + . . . , (26) 
10 



where dots denote the terms at other powers of In /3 and j3. The quantities af{fi) and 6|^(/u) depend on 
the spectrum parameters. The vacuum energy divergencies, appearing when /3 tends to zero, should 
be canceled by appropriate counterterms in the initial action of the model. This renormalization 
procedure reduces, in essence, to a replacement of the divergent in the limit /3 — )■ coefficients at 
the different structures a^(0) In" 6^(0), n < I and k < 0, hy some finite constants. These constants 



have to be fixed by certain additional normalization conditions. Keeping in mind (25), we see 
that the total renormalized average energy of the system of fermions with their vacuum at a finite 
temperature does not contain the term ag(0) in the following sense 



En 



ren. 

aO(;u) - aO(0) + (a?(/.) - a?(0)) ln(/36?(/.)) + a?(0) In MM + . . . , (27) 



where c is some constant. This property holds for bosons as well, inasmuch as the vacuum energy 
can be obtained from the high-temperature limit of £'b(/3,0) — 2Eh(2(3,0) = Ef{p,0). However, 
it does not mean that the quantity aQ(0) cannot be observed at the vanishing chemical potential. 
At low temperatures, the thermal contribution to the average energy tends exponentially to zero, 
whereas the vacuum energy is independent of the temperature and has the form ( p6|) with /i = 
and renormalized (not necessarily to zero) divergencies. 

To a ceratin extent, the above reasonings can be generalized to the multi-loop corrections with 
a regulator of the thermal form, but we leave a detailed study of this question for a future research. 



1. Expansion of the quasiclassical contribution in terms of poles 



The quasiclassical contribution to the sum (0) at g = can be also represented in the form 



analogous to (P4|). To this end, we put q = in the integral (18). Then we obtain 
^d/2-1 



T{d/2) 



dsf{s) — 2i I dslmf{s + ie) 



00 p'^-i 
dp 







2 ' 
— s 



Red E (0,2). (28) 



The integral over p reduces to the beta function. Taking the integral over s and using the relation 
(I23D, we arrive at 



00 

^2/0 = -7r'^/2^r(-(i/2)e-^^'^/2 ^S^pf - 7r'^/2-i^(l - d/2) j dss'^/^'^ lmf{-s + ie), (29) 

' \So\ 

where the summation is carried over all the poles of the function /(s). Note that this expression is 
real when the variable d is real. It is, of course, an anticipated result. In order to take advantage 
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of this formula one needs to sum the series over k, represent the result as an analytic function of 
d, and only then set d to the physical dimension. The same procedure should be applied to the 
integral over the cut. Also, if it is necessary, one can deform the integration contour to make the 
integral convergent. 

Integrating Eq. (pg|) over at that d when the resulting series in the right-hand side (RHS) 
converges, and taking into account the boundary conditions (11), we can write 



7r''/'^T{-d/2)e 



-i-nd/2 



1 



+ -7r^/'-'^{l-d/2) j ds 

\Sc\ 



s'^/2-iima;+(s) + d"is'^/2 (c^'+(s)tanh 



, (30) 



where we have integrated by parts like in Eq. (|24|). If |pfc(/i)| tends to infinity with ^ tending to 
minus infinity then the last term in the square brackets in the first line vanishes. This is the case, 
for example, when the dispersion law ijo{s) has a power-like asymptotics at \s\ — )• oo. Just as for a 
contribution from the cut to the essentially quantum part of the i7-potential, ilg, the second term 
in the square brackets disappears when the values of the function oj{s) taken on the opposite banks 
of the cut differ by a sign only and we take into account the contribution from antiparticles. 



B. Expansion effectiveness conditions 

In order to clarify the conditions of an effectiveness of the expansion ( p4[ ) , it is useful to provide 
it by a simple visual interpretation. As we have already mentioned, every term of the series over 
poles in Eq. (p^) substituted to (|6|) can be thought of as the Euclidean propagator of a particle 
with the mass ^epk and the rate of decay Imp^ in {d + 2)-dimensional space. The sum entering 
the expansion (|2^ ) can be interpreted as a sum of wave-functions (taken in the origin) of these 
"particles" propagating from the lattice points with coordinates and the initial phases 
This lattice is embedded in a (ci -|- 2)-dimensional Euclidean space, whereas the Gram matrix of the 
lattice spacings specifies the metric gij. 

Now it is evident that the sum over poles converges rapidly if the lattice spacings are large 
enough, i.e., the elements of the Gram matrix g^j are large, and the "screening parameter" Imp^ 
tends rapidly to infinity with the number k. Then, with a neglible error, we can break off the series 
in k and q retaining only several first terms in them. We see from Eq. ([l5| ) defining the singular 
points that the imaginary part of pk grows with k for oo{s) ~ s", a > 0, at large s, it growing 
slower at the greater a (with the exception of a = l/4n). If a;(s) ~ exp(as), the imaginary part 
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of pk decreases. Thus, if the dispersion law a;(s) grows with s more rapidly than any power of s 
then the expansion (24) cannot be used to make an accurate estimation of the sum (DQ. Further, 



we shall assume that Imp/^ and \pi^\ increase monotonically when an absolute value of the number 
k increases. 

To derive a rough estimate of the expansion effectiveness conditions, we require that a term of 
the series in k with the number k^ will be much greater than the next term of this series starting 
with some small number k^. Let us use the series representation of the Hankel function at a large 



argument j45 



p) « ^-^e-^-^p-W2-m T r M" ^^'"''^' r(. + n + l/2) 



„_o / g'^+^+^Z^ nW{u - n + 1/2) 
This is an asymptotic expansion. If we take only the first M terms of this series, we make an error 
with an absolute value lesser than an absolute value of the first discarded term. The reminder 
term has the same form as a general term of the series taken at n = M + 1 and multiplied by a 
number with a modulus lesser than one. For a half-integer z^, this series terminates giving rise to a 
closed expression for the Hankel function in terms of elementary functions. Using this asymptotic 



expansion, we have from (24) the condition of a rapid convergence 

> 1, (32) 



. d 

2Trqm lm{pko+i - Pko) - -^n 



Pko+l 



where := mm{q), q ^ 0. This condition constrains the values of the metric components gij 
which characterizes the sizes of a system. In particular, this condition implies lirq^^i^Pko+i > 1) 
what allows us to break off the series in q for the contributions from the poles with numbers ko + l 
and higher. In order to terminate the sum over q for the contributions from the poles pk, \k\ < ko, 
it is necessary to demand, in addition, 

2TTqmln\pk > 1, for \k\ < ko. (33) 

We shall show below how to obtain an expansion of the logarithm of the partition function in that 
case when the last condition violates for a few number of poles. Notice also that the expansion 
effectiveness conditions can be used to specify a domain of values of the system parameters in which 
the oscillating non-extensive corrections described by the pole terms appreciably contribute to the 
0-potential. This domain is at a saturation border of the inequality (|33|). 



^ If the function uj{s) grows rapidly with s, a good approximation for the partition function can be obtained directly 
from its definition (j^) discarding all the higher terms of the sum. Analogously, we can obtain a good approximation 
for the partition function in the case of small metric elements gij. 
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As one can see from (31) and the interrelation between the Macdonald and Hankel functions 



(|2lD , the sum over q in the contribution from the cut to the expansion (|24|) can be broken off if 

2^<7m|sc|'/' > 1. (34) 

The method of obtaining of a rapidly convergent expansion in the case when the branch point Sc 
approaches zero will be considered below. 

The expansion effectiveness conditions can be written more explicitly if we know a concrete form 
of the dispersion law. For example, in the relativistic case, the poles defining the expansion of the 
il-potential becomes 



Pk 



\/ {fi + iu)kY ~ "T-^) Impfe > 0, k £ Z. (35) 



Recalling our choice for the branch of a square root entering the dispersion law, we see that /u should 
be positive. When we decrease the value of ^, the poles in the s-plane approach the cut (—oo, —rn^] 
from the above and below in a symmetric way. At /i = all the poles lie on the cut and have to 
be taken with the factor 1/2 in the expansion (p^). At /i < all the poles leave the "physical" 



sheet and the sum over them in (24) vanishes. From the physical point of view, this corresponds 
to the fact that, at negative chemical potentials, particles and antiparticles switch their roles - the 
probability for a particle to appear in the system becomes exponentially suppressed. 



At sufficiently small temperatures (m and are much greater than uiko), the condition (32) 
approximately reads as 

(36) 

2'K'^qm{uJko+i + Wfcjm^ > (m^ - fi'^f''^, < m. 

As for bosons, we need to use the second condition. In the Boltzmann limit, fi < m for fermions 
as well (see, e.g., pg|]), although, for a degenerate Fermi gas, the nonrelativistic chemical potential 
/^nr = /i — m is positive. At high temperatures (m and fi are small in comparison with oJko), the 



condition ( p2D in the leading order takes the form 



^TT^qm > 1 + ^ In 



(37) 



This condition is, in fact, the expansion effectiveness condition (24) for the oscillating part of the 
partition function in the case of a massless quantum field at zero chemical potential. 
For the nonrelativistic dispersion law, the poles look like 



Pk = y/ji + iujk, Impk > 0, k eZ. (38) 
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At low temperatures (/i is much greater than UJko), the condition (32) in the leading order reduces 
to 



/" > 1; 



vr^gm(t^fco+i + ^fco) ^ 2|^|^/2, /i < -1. 



(39) 



At high temperatures, the condition (32) becomes 

1 + hi\ujkQ+i - ^fcol"'^^ A 



(40) 



where we have taken into account the first correction in jjL. One can also derive these conditions 
(without the first correction) as a nonrelativistic limit of ( |36[) and 



Resummations 



1. Summation over poles 



Let there exist a sequence Sa, ol = N,cc, among the poles such that the distance between 
adjacent members of this sequence \pa+i — Pa\ tends to zero with an increase of the number a, 
whereas Im pa tends to infinity. Then, the estimate of the ri-potential, which is obtained by taking 
into account only the first several terms of the expansion (24), can be improved summing the 
sequence Sa by the Euler-Maclaurin formula 

where the integration is along the contour which connects the poles Sa starting from the A^-th pole. 
In consequence of the symmetry (p!^), the sequence of conjugate poles s* also contributes to the 
expansion (p4|). Keeping only the leading integral term, we obtain for a contribution of these two 
sequences 

- Vff Re r -u;'is + m')dsh'~^l^^{q, s^'^) + ... (42) 

The integral in the Euler-Maclaurin formula can be taken for a nonrelativistic dispersion law when 
uj'{s) = 1 and the distance between the poles tends to zero at A; — )■ oo. This holds for any 



polynomial spectrum as well. Formula (42) resembles the contribution from the cut entering the 
expansion (j2^). It is not surprising since the poles merge and constitute a cut when a distance 
between them is small. 
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2. Summation over q in one dimension 



At the large argument, 

iTrq\p\ > 1, (43) 

the Hankel function admits the asymptotic expansion (^) . If the summation over in the series 
reduces effectively to a one-dimensional sum, viz. 

q = 9n^\q^\=qm\q^\, (44) 

then we can sum the series over q 

X [Li,+„+i/2(e2-('''"^'+"^))+Li,+„+i/2(e2-('""f-'^i))] , (45) 
where we have introduced the polylogarithm function 



The above described situation occurs, for example, in a one-dimensional problem or in the case 
of a thin film. In the latter case, the summation is carried over that quantum number which 
corresponds to the small dimension of a film. The contributions of transverse dimensions to the 
oscillating part of the O-potential are strongly exponentially suppressed. The obtained expansion 



can be also used to evaluate the contribution to the partition function of the integral (42) 
coming from the Euler-Maclaurin formula or the integral representing the contribution from the 
cut. In order to apply it in these cases, the argument of the Hankel function should be large enough 



on the integration curve, that is, either the condition (43) is fulfilled for the pole pi\[ at q = q^n or 
the Hankel function has a half-integer index. 

By the same way, we can sum over q along any vector on the lattice in the (/-space. This allows us 
to rewrite the whole sum over q as the sum over all irreducible vectors, whose components have not 
a common divisor excepting unity, in the g-space. As a result, we obtain a sum of the expressions 
of the form ( ^5|) over all the irreducible vectors of an increasing length. This resummation is an 



effective tool in the "one-dimensional" case or in the case when the condition ( 33 ) is weakly violated 
for the given pole p^. 
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3. Contribution of a pole approaching the real axis 



Let us now study the case when the condition ( |33D is not satisfied by some pole and the series in 
q cannot be terminated. This case takes place when that pole approaches the real positive semi-axis 
in the s-plane 

SQ = h + ie, b>0, (47) 

and e is a tiny quantity. Besides, we can assume that is not small since, in the opposite case, 
the partition function in its initial form (|5|) rapidly converges and can be simply estimated by the 
first several terms of the sum. Usually, the poles approach the real axis at low temperatures (see 
Eq. ([l5|)) or, for the Bose-Einstein distribution, at a small chemical potential like in the case of 
a condensation. In spite of the fact that, in this section, we are about to investigate the pole 
contribution, all the analysis presented below is also applicable to the Macdonald function entering 
the integral over the cut in (^), when the branch point Sc is in close proximity to the origin. 

Consider separately the contribution from the pole ( ^7|) to the oscillating part of the partition 
function (^ 



^J^e'^''^^^'h'jl{q,Po). (48) 



-ITT 

1 



Recalling the integral representation of the Hankel function |48|, 

idt(ii)-i-^/2gi(-2gVt+«ot)^ Ijnsg > 0, Red < 0, (49) 

the recurrence relation ([To|) and the symmetry property (^), we see that h^\q,p) is an entire 
function of the variable z/ at g 7^ and p M. As long as the series ( p8| ) converges uniformly with 
respect to d G C, it defines an entire function of the variable d. It is convenient to consider this 
series at Re d < when we can add the term g = to the sum. In this case. 



c := Cos - ir^'^M-Sof'^n-dm = -zn^^e'^^-^"' h'jl^{q,po) = - ^ vr"^r(M)(p2 - so)-^ 



(50) 

where Pi := pi + aj and fi := {d — d)/2, while d is the physical space dimension. Regarding fj, as 
an arbitrary parameter, we shall derive the rapidly convergent expansion of the contribution c(/u) 
from the one pole. Then we shall set fi to zero. For a definiteness, we shall assume that e > 0. The 
expansion for the negative e can be obtained by a complex conjugation. So, 

rX poo 

c=-7r'^/2^V/ idt(it)-i-'^/2gi(.2gVw+2™,,0_^-MV / idt{itr-^e-''^'''-''\ (51) 

g Jo p JX 
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where A is some complex number. It is useful to choose A to lie on the imaginary axis, A = — i|A|, 
although it is not obligatory. In the first integral, we pick out the term q = 0, 

c? := -7r''/^^{-sof'[n-d/2) - r{-d/2, -\X\so)], (52) 



where r{u,z) is an incomplete gamma function |48| 



r{u,z) = T{iy)-y^ -. (53) 

The first term in the expression for is exactly the contribution to the sum ( |50D at g = 0. This 
term does not enter Cqs- The contribution from q ^ to the first integral in Eq. ( [5l| ) can be written 

as 



^y^^ \s'(t)dt) s'it) 



|A|'^/2 

' ' q n=0 



, (54) 



t=i 

where S'{t) := 7r'^q'^\X\~^t~'^ + sqI'^I- The expansion in the first line should be exploited when 
the quantity 7r^g^|so| is small, while the asymptotic WKB-expansion in the second line is useful 
when |5"(1)| is large. The second integral in Eq. (^Tj) is expressed through the incomplete gamma 
function 

Coo := -TT'^' Y.^P' - sor^r (^, |A|(P2 _ so)) . (55) 
p 

Let us remind some properties of the incomplete gamma function which are needed for our 
subsequent investigation. The incomplete gamma function T{i',z) is an entire function of at 
z ^ 0. In the z-plane, we choose its branch with a cut along the negative real semi-axis for N. 
The incomplete gamma function T{i',z) is singular in the point z = at Keiy < 0. Its absolute 
value grows exponentially at Hez — )• — oo and exponentially goes down to zero at Re 2; — )• 00. In 
particular, we see that the contribution from the one pole to the oscillating part of the sum (P), 

Cos = ci + Coo + c? + 7r'^/'V5(-so)''/'r(-d/2), (56) 

is an entire function of v as expected. Now, we can put fi = 0. 

The parameter A is an arbitrary number and resembles the massive parameter fj, of the dimen- 
sional regularization of quantum field theory. We choose the value of the parameter A so that the 
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sums over q and p in the expressions for ci and Cqo can be broken off with a neghble error. In order 
to terminate the sum over q in (|5^), it is necessary to demand 



> 1- (57) 

Then the term ci is exponentially small. At the same time, A should not be very small since, in that 
case, the sum ( |55[ ) over p converges slowly. The main contribution to the sum over p comes from 
the points in a vicinity of the "resonance" = b and from the exponentially growing contributions 
at a large negative real part of the argument of the incomplete gamma function. By use of the 



expansion (53), it is not difficult to estimate a half- width of the maximum of real part of the 
incomplete gamma function entering C2, as a function of P^, at fi = 0, 



A = 2^p^.--'^ (58) 

where 7e is the Euler constant. We are interested in the real part of the contribution because the 
poles appear by complex conjugate pairs in the expansion (^) such that the imaginary part is 
canceled out. To chop off the exponentially growing terms (at small P^), we additionally assume 

|A| < x/b, X « 0.4. (59) 
The last number is taken directly from the plot of the incomplete gamma function at = 0. If the 



conditions (57) and ( p9| ) are fulfilled then we have to retain only the first several terms in the sum 
over q which are near to the point q = 0. As for the sum over p, we should keep only such terms 
that 

p2 e (6- A/2,6 + A/2). (60) 

If the metric gij is approximately isotropic (all of its eigenvalues are of the order of g^) and e is 
small then the expression for Cqs contains only a few number of relevant terms. 

The case of a strongly anisotropic metric needs a more cumbersome procedure. In this case, 
an immediate application of the above formulas is not very effective since there are a lot of terms 
which satisfy the condition (|60|). These are the terms that correspond to the vectors Pi differing by 
the components along the large principal axes of the metric gij. This problem is also easy to see 
from the expression ( p5[ ) when |A| is of the order of the least eigenvalue of the metric. 

Now we show how to modify the above procedure to the case of an isotropic metric. Let the 
eigenvalues of the metric gij constitute a hierarchy 

Ai < A2 < • • • < An. (61) 
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For brevity, we shall suppose that the metric is diagonal. Then, we part the integral (^Tj) according 
to this hierarchy 



c = — vr 



r /"Ai /■A2 roo 

''/^^V / idt+ idt + ...+ idt 

„ \_J0 JXl JXn 



i(7r2qr2/t+*ot+27raig*) 



(62) 



where Aj are complex numbers with moduli being of the order of the respective eigenvalues Aj. As 
before, we assume that Aj = — i|Aj|. We have already known the first and the last integrals in 
It remains to consider a contribution of the A;-th integral 



It is useful to denote by all the indices q corresponding to the subspace of the metric eigenvectors 
with the eigenvalues Aj < A^-i, and by q^ all the rest indices q. We also denote by Rk a number 
of indices q"" and assign by definition Ri := 0. Then 



ga 



(d-Rk)/2 9_y^ ^^^_i_(d-Rk)/2^-7T^qyt-{PS-so)t-27riaAq^^ (g4) 



where ga := det gab, q\ '■= dABQ^Q^ , and so on. Taking into account the range of variable t, it 
is now easy to see that the sum in the last formula can be terminated with a neglible error. A 
contribution of the zero mode has the form 



cO=_^id-R,)/2l9_ 



\(d-Rfc)/2 



[T {{Rk - d)/2, \Xk-i\{P! - so)) - r {{Rk - d)/2, |A,|(P2 _ so))] . (65) 



This expression is a monotonically decreasing function of P^. With an increase of the indices Pa, it 
tends exponentially to zero. Hence, we can retain only the first leading terms in this sum for which 
P2 « 0. 

An asymptotic expansion of the integral at q^ can be obtained by the WKB-method (see, 
e.g., |[50|] ). This method leads to the five cases corresponding to different values of the indices q^ 
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and Pa- They are 



Cfc 



ga 



-2'rriajiq 



(vr/iA.i)('^-«'=)/2E:r= 



1 t-l-(''-«fc)/2 



S[(t) dt 



=Si(l) 



f{d-Rk)/2{qA, \/Pl - So), 



f{Rk-d)/2 ( V^'a - So, g^) , 

-(^/iA,_ii)('^-«^-)/2E:r=( 



-1 d 



^-l-(d-Hj.)/2 



vrM/|Afc| + |Afc|(6-P2)>l; 
vrM/|Afc| + |Afc|(6-p2)«0; 



t=i 



7r2<;2/|Afe_i| + |Afc_i|(6-p2)«0; 
e^^W, vr2g2/|Afc_i| + |A,_i|(6-p2)<_i 



where ky{q,p) is the Macdonald function defined in Eq. (^). We also introduce the notation 

= -k4 + l^'^K^o - ^')*' ^2(t) = + lAfc^iKso - Pl% 

W \\k-i\t 

ri/p 



JO 

The function fuiQ,p) is not a cyhnder function and apparently cannot be expressed in a simple 
form (see the remark on p. 313 in ||4^ ). One can think of this function as a "half" of the Macdonald 
function in the sense that 

Mq,p) + f^,ip,q) = 2K{q,p). (68) 



At large arguments, the function fu{q,p) possesses the asymptotic expansion analogous to ( |3l| ) 

f(an)-^(PY e-^-'^P 2"(-l)»-feC^ r(l - z.)r((n - k + l)/2)r((n + k + l)/2 + s) 

JAQ,P)-y^J Z^Z_>Z^ (47rgp)(n+fc+i)/2+s nisir{l-u-n)r{{n-k + l)/2-s) 

71^ — s — Aj — 

(69) 

In increasing with fixed q"^, the above mentioned cases are realized successively from top to 
down. It is clear from the representation of the integral ( |64|) that its value rapidly decreases when 
P^ increases. One can see this from the explicit expressions (p^) as well. In other words, to obtain 
an estimation of the contribution of this integral to the il-potential, it will be sufficient to choose 
the first case (at the given value of b) from the indicated ones and take a few first values of the 
indices q^ and p^. 

In summary, the contribution to the oscillating part of the partition function from a pole ap- 
proaching the real axis reads as 

n 

Cos = ^^(q + c°) + Coo + vr'^/2^(-so)'^/2r(-d/2), (70) 
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where we should take A to be equal to A„ in the expression for Cqo- The value of A„ needs to meet 
the condition (|59|). Otherwise, in the sum (|55|), we have to take into account the contributions not 
only from the indices pi satisfying (|60D, but also from that indices pi for which is small. 

D. Quasiclassical contribution 

1. Nonrelativistic limit 

Let us turn to the zero mode g = of the integral (quasiclassical contribution) . On integrating 
over angles, it can be cast into the form 

= r + e---(-^+™^)). (71) 

Now we introduce a useful notation 

a;(p^ + m^) - Wo =: y(s), ojq := ujim^), uio := ujq - fi, (72) 

where s := p^. In order to get rid of the logarithm in the integral ([n|), we integrate it by parts and 
substitute a series representation of s"(y) to it, 



oo 
fc=0 



1 d^^ 



y'{s) ds 



y{s) 



(73) 



Making use of the definition of the polylogarithm function (^6|) , we arrive at the expansion (for the 
relativistic dispersion law see, e.g., [ p7[| ) 

The analogous expansion for the Bose-Einstein distribution is recovered by a change of the overall 
sign and the sign in the argument of the polylogarithm. Since y is proportional to /3, the expansion 
( |7^ is carried out in inverse powers of j3 and provides the low-temperature asymptotics of the 
quasiclassical contribution to the ri-potential. In particular, it terminates on the first term for the 
nonrelativistic dispersion law. 

2. High-temperature expansion 

In order to obtain the high-temperature expansion, we shall proceed in the analogous way. On 
integrating by parts in the integral (ffll), we have for different statistics 
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where we have shown exphcitly a dependence on /3 and supposed that uj{s) is a monotonically 
increasing function. Our aim is to expand this integral in increasing powers of /?. First, we expand 
a numerator of fraction entering the integrand. Assume that, at sufficiently large \y\, the following 
representation holds 



(—^\p 



>0 p! 

p=0 ^ 



oj'{s) ds 



-\ a 

S 



(76) 



where 7 > and uj ^{y) is an analytic function of 7. Then, sequentially expanding the expression 
in series, we obtain 

(,r'(B-^„^„^m^\'^/'- V (-1)' 4^'-'rid/2 + i)r(7(d/2 -k)-p + i) 

[u [p y + i^) m) -^2^^ r{d/2 + l-k)T{j{d/2-k)-l-p+l) 

X m^''^\l3~^yy'^'^''^-^'^-^-P. (77) 

This power series converges at \y\ > /3i? at sufficiently large R independent of /3. 
Further, we need to use the well-known expansion of the incomplete ("-function 



/ 

J a 



n=-l 



^^^4-7 = (1 - 2'"'^)r(^)C(^) - E(i - 2 



i+„,(-irc(-nK 



(78) 



+ l ^0 r(n+l)(i. + n) 



In the bosonic case, this expansion is valid at \a\ < 27r, while, in the fermionic case, it holds at 
\a\ < TT. The functions in the RHS of the equalities are entire functions of v. Let us part the integral 
( |75D into two: [/3tDo,/3-R] and [/3R, +00); and denote these integrals as ii and 12, respectively. Then, 
at sufficiently small /3, we can apply formulas (|78[) to the second integral ^2- Substituting the 



expansion (77) into 12, we integrate the series term by term. As a result, the whole expansion of 
the integral 12 splits naturally into two contributions from the first and second terms in Eqs. (|78D, 
respectively. In the first contribution to 12, the expansion is carried out in increasing powers of /3, 



there being a finite number of terms at any fixed power of /3. Indeed, in the bosonic case 

kW.T{d/2 + l-k) 



i2 = ^"^v~9 E m\Tui2 + i-k\ nim-k)-p+mi{d/2-k)-i-p+i)+^m), 

(79) 

where {p{f3R) denotes the second contribution. We also explicitly marks its dependence on the lower 
integration limit in ([f5[). This last term has to be resummed. From Eq. ( [7^ ) it is not difficult to 
see that every term of the series in n in the expression for ip{f3R) can be written in the integral 
form 
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T{d/2 + l) r(n + l) 



where d and 7 are chosen so that the integral converges and then continued by analyticity to their 
initial values. Bearing in mind the Taylor expansion of the function (e^ — in y in a vicinity of 
zero, the expression for the integral ii can be cast into the form (pO|). Matching the parameters d 
and 7 in the both integrals, we can combine them into one. This gives 



^ r(n + i) ^ r(fi/2 + i) ^ ^ 



n=-l 



For the Fermi-Dirac distribution, the series (|79D and (|l]) have the same form as in the bosonic 
case. The difference consists in the overall sign, powers of two typical for the fermionic distribution 



and other summation limits: n = 0,oo (cf. Eqs. ([78|)). Thus, in both cases, the problem reduces 
to a finding of the integral 

-^'-2im^,) = l^dyiy-,r ^ r(d/2 + 1) =1 rWTI) " + ^-(^ + " , 

(82) 

in the form of an analytic function of the complex variables d and 7. It is clear that this integral 
converges at 

Red > -2, Re{jd) < -2n - 2. (83) 
The factor m'^ is introduced to Eq. (|82[) just for convenience. 



If the parameters m and ^ are arbitrary, the expression (82) for o"^ does not contain any small 
parameter and cannot be evaluated or expanded. Therefore, we need to specify a form of the 
dispersion law uj{s) in order to proceed to a further investigation. For example, for a homogeneous 
dispersion law, 

Lo = s^'\ X- = (84) 



the sum over p in the expression (79) stops at the first term and the integral ( |82D reduces to 



^ ^2(n+l)/7 /"°° _ z^^^^L^^^ 

Ji r(d/2 + i) 



m 



2( 



V -n- l)/7 - d/2) T{n + l){-zY 
^ ^ r(fs-n-l)/7 + l) sT(n-s + iy ^ ' 



((s-n-l)/7 + l) sT{n-s + l) 



=0 

where z = fim~'^/"' . This series representation is valid at \z\ < 1 and also at the nonnegative integer 
n when the series terminates. The high-temperature expansion for the Fermi-Dirac distribution 
depends on a'2 with nonnegative integer n only. So, in the fermionic case, the high-temperature 
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expansion of the integral ( [75D reads as 

I', = -'"^aY: E(i - 2'""^'^'"'^)^iir^^|^:T3fe^r(7(d/^ -k) + m^{d/2 -k)-i + i) 



1=0 



kj-k- 

'■ (86) 



Note that a "half" of the contributions from the second term in the square brackets vanish due to 
the property of the C-function. These are the terms at the even positive I. 

In order to obtain the high-temperature expansion for the Bose-Einstein distribution law, we 
need to evaluate the integral 

= A r(ci/2 + i) =^ g r(//, + i) ^''^ 

This integral is an analytic function of z with a branch cut discontinuity along the part of the real 
axis z > 1. The discontinuity on the cut is equal to 

aj\z + ie) - a^\z - ie) = 2^. ^^^^^^^^ . (88) 

The above series representation holds for \z\ < 1. It allows us to obtain the analytic continuation 
of this integral to arbitrary complex values of d and 7. If 7 is a natural number then this integral 
can be written as a finite sum of the Gauss liypergeometric functions 

7-1 



n=0 



At 7 = 2 (relativistic dispersion law), we have different equivalent representations 



= rfi^irf/2) ^(l^ 1 - (1 + z)/2) = -^i^L-^F{l/2, -d/2; 1 - d/2; 1 - z 



2\ 



^r{-d/2){l-zY/^+T ((1 - d)/2) -^F{1, {l-d)/2; 3/2; z') = 2''/'T{-d){l-zY^^P^J/l{-z), 



(90) 



where Puiz) is an associated Legendre function of the first kind |51| single- valued and regular at 



\z\ < 1. For an arbitrary 7, it is easy to obtain the following recurrence relation 

^,-.(./2+i)-v^-i(^) = 9.M^/2+%,-i(z)]. (91) 

To derive the expansion of (^^^{z) at \z\ > 1, we write it in the equivalent form 

da J (irE(x'^-l)'^/2g-a(x-^)^ Rez<0. (92) 
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Now we expand the preexponential factor in x and integrate the series term by term. Then, 
by the use of the expansion for the incomplete gamma function (|53|), we expand the obtained 
preexponential factor in a and again integrate the series term by term. Thereby, we come to 



(-1)'= 



fe=0 



klT{d/2 + l-k) 



7r(-z)7W2-fc) 
sin7r7(A; — d/2) 



ra=0 



-1-n 



7(d/2 -k)+n + l 



(93) 



The series in k in the last term in the square brackets is easily summed and expressed in terms of 
the beta function. As a result, we have 

(-1)' 



k=0 



kT{d/2 + l-k) siuTTjik - d/2) 



r(i -{k + 1)/7) 



(94) 



If 7 is a natural number, we deduce the representation 



sm{TT-fd/2) T{d/2 + l) ^ ^ V(\-kh) 

k=l 



' t ^^~t'^~:!'^ ^''ni. k/r, kh + d/2 + 1; .-). 



(95) 



Thus the high-temperature expansion in the bosonic case becomes 
°o r / ^2\k,Am{k-d/2)+l 

n = -'"^aY: E . ni{d/2 -k) + mj{d/2 -k)-i+i)+ 



1=0 



k=0 



kUT{d/2 + l-k) 



+ 7- V'^(/3m2/T)'+^C(-0 E 



k=0 



(-l)^z'-^T(-rf/2 - (fc + l)/7) 
k\{l-ky. r(l - (A; + l)/7) 



(96) 



It is clear that the obtained high-temperature expansions (^6|) and (^) are straightforwardly gener- 
alized to the dispersion laws in the form of a finite sum of homogeneous functions. In the expression 
(|96|), as well as in the high-temperature expansion for the Fermi-Dirac distribution ([8^), we can set 
d to be equal to the physical dimension and, henceforward, we shall denote the physical dimension 
as d. Then the two terms in the brackets and the last term possess the poles in the 7-plane which 
are mutually canceled out. 

In the fermionic case, the singularities arise when the physical degree of homogeneity 7 of the 
dispersion law is a rational number of the form 

2p + 2 

lor OQQ a; 

(97) 



7 
7 



2q + r 
p+l 



+ r 



for odd d; 
for even d; 



where p and q are nonnegative integers. Denoting by li and ki the indices in the first term of the 
high-temperature expansion ( p6| ) and by I2 and /c2 the indices in the second term, we see that the 
cancelation of singularities occurs between the terms with 



li — I2 — k2, k2 < I2', 



k,--r 



d k2 + l 



7 



(98) 
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d 



+ q + n{2q + + n{2p + 2] 



n = 0, oo, 



ikuk2)= { -+n{q + l),n{p+l)-l] , 



n = 1, oo, 



for odd d: 



for even d. 



(99) 



In the bosonic case, there are extra singularities coming from the pole of the ^-function in unity. 
They are canceled out by the singularities of the last term in (|9^ ) . Expanding this last term in the 
series ( p7| ) and denoting the corresponding summation index as I3, it is not difficult to understand 
that the cancelation occurs when 



h = h; li=^id/2-h), ki = 0,[d/2], 



(100) 



that is, when 7 is of the form (|97| ) and 
'd-1 



{ki,h) 
{ki,h) 



q-n{2q+l),{2n + l){p+l)] , n=0, 



Aq + 2 2 



for odd d\ 



d 



n{q + l),n{p+l) , 



n =0, 



d 



2(g + l) 



for even d. 

(101) 

Note that the terms proportional to ln(/3m^/'^) may appear in the high-temperature expansion for 
homogenous dispersion laws only with the degree of homogeneity of the form (^). Since only in 
this case we need to expand the powers of /3 and m in the expressions ( p6[ ) or ( p6[ ) into series in 7. 

As an example of the cancelation of singularities in (p6|), we consider the bosonic high- 
temperature expansion with the natural degree of homogeneity 7 G N. Besides, we are only 
interested in the singularities of the type ( |100| ) when, upon the cancelation of singularities, we 
should sum an infinite series in z to obtain a closed-form expression. Then 

= ^"'^9 E ^-iSMi^—^^nim -k) + i)ai{d/2 -k)-i + i)+ 



l.k=0 



kW.T{d/2 + l-k) 



1=0 
ffd/2~l 



' r(^/7 - d /2) 

1 r(s/7 + 



A:=0 

7 



(-l)^z'-^- T{-d/2 - (fc + l)/7) 

k\{i-ky. r(i - (ik + i)/7) 



-^"/^Vff E kiT{d/2+ 1 - k) [^^/^"^'^'-(^E + m/2 - k)))+T'mi2+i-k)-^i:{k+i)) 

(102) 

where ^/>(x) = r'(j;)/r(x), the primed sum over I denotes a summation over those I at which the 
argument of the (^-function is not unity, whereas the prime at the sum over s says that we throw 
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away the singular terms arising from a nonpositive integer value of the argument of the gamma 
function. In the odd-dimensional space, the above high-temperature expansion holds only for even 
7 where there exist singularities of the type (100). Specific values of the Gauss hypergeometric 



function can be found, for example, in |52|. In the case of relativistic and nonrelativistic spectra, 
all the hypergeometric functions entering the last formula are expressed in terms of elementary 
functions. 

In conclusion of this section, we write out the linear in /3 terms of the high-temperature expansion 
for the fermionic f]-potential at a vanishing chemical potential. As we know, these terms are of 
a particular interest since they have the same form as a finite part of the one-loop contribution 
to the vacuum energy (p5[). Starting from the general formula it is not difficult to find this 



contribution. Tending 7 to its physical value 7, we have 



;; = ,./.^^„«/,IIz|^ + ... (103) 



when d/2 + I/7 is not natural. 



2r((i/2-h 1/7 + 1) ^ ' 

when d is even, but 7^^ is natural, and 

^'^ ~ 2r(i - i/7)r((i/2 + 1/7 + 1) 

In + 7E + 7"' [2V'(2 - 2/7) - ^(3/2 - I/7) - ^(d/2 + I/7 + 1)] U . . . (105) 

when d is odd and 7"^ is half-integer. The last formula gives, in particular, a finite part of the 
quasiclassical contribution to the vacuum energy for particles with the relativistic dispersion law in 
a three-dimensional space. 



III. EXAMPLES 

A. Massless particles 

1. Zero chemical potential 

In this section, in order to illustrate the main features of the developed general procedure, we 
consider a simple model of a gas of free massless particles confined to the rectangular parallelepiped 
with sizes {Lx, Ly, L^) both at zero and non-zero chemical potential. We start with the case of a 
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vanishing chemical potential Such a model describes, for instance, a gas of photons in an ideal 
metal rectangular box. Notwithstanding this model is profoundly investigated (see [^, |l^ and some 
recent papers [pC|-p3|, ^-^|), its consideration will be instructive because of its high degeneracy. 



We shall derive new rapidly convergent expansions for the rj-potential applying the resummation 



formulas found in Sec. lIC. For the sake of conciseness, we shall study the fields subjected to 
the periodic and Dirichlet boundary conditions. For the model at issue, the partition function 
corresponding to the quantum fields with the Neumann boundary conditions can be reduced to the 



partition function of the fields with the periodic boundary conditions (see, e.g., |21, |2^). 
The "metric" is defined in this case as 

gij = l^'^diag{L'^, Ly, lT = Trl3ch, (106) 

where It is a thermal wave-length of a massless particle. This quantity is a characteristic scale of 
the model. At the temperature /?^^ = 1 K, it approximately equals to 0.72 cm. The roots of Eq. 
( |T5|) take a simple form 

Pk = fJ- + iuJk, pI = -ji + iuJk, k>0. (107) 



As we have already noted, the poles (107) lie on the "physical" sheet of the function u}{s) only 
if /X > 0. At the vanishing chemical potential, they contribute to the oscillating part of the Q- 
potential with the factor 1/2, the integral over the cut being understood in the sense of principal 
value. Further, we shall assume that /x = 0, but, in some formulas, it will be useful to put /i to zero 
only in the finial result. 

Let us first suppose that a summation over the quantum numbers pi (do not confuse with the 
roots pk of Eq. ( |107| )) is carried over the infinite ranges (— oo, c«). Then the oscillating contribution 



to the logarithm of the partition function in d dimensions reads as 

oo 

±(3nf = -m^Y^Y.^^^/2i<l^iuju), (108) 

q k=0 

where minus corresponds to bosons and plus is to be taken for fermions. In the bosonic case, the 
term at A: = comes with the factor 1/2. As long as Repfc = 0, the oscillating contribution does 
not, in fact, oscillate since the Hankel functions turn into the Macdonald ones with real arguments. 
However, we shall call it oscillating to distinguish from the other contributions to the ri-potential. 
By the use of the asymptotic expansion of the Hankel function (|3l|), we can sum the series in k. In 
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the fermionic case, it becomes 



f27T\ ^ (87r2g)-T ((d + l)/2 + s) 



27r^\qJ s\T{{d+l)/2- s) 



X [Li,_(,_i)/2 (e-^-''?) - 2-('^-i)/2 Li,_(,_i)/2 {e-'-"")] ■ (109) 



If the condition (37) is satisfied then the series in q can be broken off taking into account the 
first several terms only. For odd d, the series in s terminates, all the polylogarithms entering the 
formula are expressed in terms of elementary functions, and we obtain the exact expression for the 
il-potential. As for even d, the series in s can be also terminated with a neglible error. The optimal 
number s, at which the series should be broken off, is the number of the series term with a minimal 
absolute value (see the remark after Eq. (^I])). The oscillating contribution becomes significant 
only when the characteristic scale of the system is of the order of It or smaller. On the scale It, a 
relative contribution of the oscillating term grows exponentially with decreasing sizes of the system. 
In accordance with the general formula (24), the contribution to the il-potential from the cut is 



written as 



Q 

'-d 



^ E' =■ ((. + 1)/2) 1), (110) 



where Ce(5>'^) ^^^^ homogenous Epstein ("-function. The series defining this function converges 
at those values of and d which we need. However, it is a slowly convergent series. One can resum 
this series to write it in the rapidly convergent form. This can be done, for instance, by means 



of the nice formula from |18, |2^, but here we apply the method developed in Sec. II C 3 . If 



the eigenvalues of the metric constitute the hierarchy ( |61| ) then, in the notation of Sec. II C 2 , the 
resummed homogenous Epstein (-function becomes 



n+l , .,y-{d-Rk)/2 _ .u-{d-Rk)/2 



p k=l ^ 



u-{d-Rk)/2 

+ Y^pt'''-^'' [^{^ - {d- Rk)/2, Xk-ipl) - Tiu - {d-Rk)/2, X^pD] 



Pa 



+ 5^'(^(Za)"+^* [r((J-i?,)/2-z.,7rM/Afe)-r((J-i?fc)/2-z.,7rM/Afc-i)]|, (HI) 

where Afc are arbitrary real numbers of the order of the eigenvalues A^, and we put by definition 
Aq = ^n+i '■— 0) ^CK' ^^id 9o ■= 1) Rn+i = d. An approximate equality means that we neglect 
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strongly exponentially suppressed terms. In the massless case, the "intermediate" integrals (64) fall 
into the third case of Eq. (H). Recalhng (|l]), it is not difficult to see that the argument of the 
Macdonald function is large and therefore it gives a neglible contribution. It is those terms which we 



casted out in (|111|). The contributions in the second and third lines of Eq. (Ill) are exponentially 



suppressed. Of course, in view of (pl|), the expansion (111) admits a further simplification, but we 
leave it intact. 

In considering massless particles at the vanishing chemical potential, we restrict ourself to the 
case of an approximately isotropic metric gij (all its eigenvalues are of the same order). Therefore, 
we have the exact expression 

~ ^ T{u/2)q- ^^TY WJW^' r(z./2 + l)+ ^ {u-d)T{u/2y 

(112) 

where A is an arbitrary parameter. It is useful to take its value be of the order of Then, 
making a neglible error, the sums in the RHS of (112) can be broken off retaining the first several 



terms only. Notice also that if and Ly are much larger than and It then the series in 



entering (110) becomes effectively one-dimensional over the vectors directed along the z-axis and 
is expressed through the Riemann (^-function. In this case, the average energy following from ( |110D 
coincides with the well-known result for the Casimir energy of a one scalar degree of freedom: we 
must also divide the expression ( |110D on 2'^ to take into account the Dirichlet boundary conditions. 
An expression for the quasiclassical contribution to the ri-potential can be easily deduced from 



its integral form. Nevertheless, we shall derive it from the representation (30) as a sum over poles 



in order to demonstrate the method. The contribution from the cut in formula (^0|) is zero for a 



massless relativistic dispersion law. An analog of the Casimir term in Eq. (30) is proportional to 
m'^^^ , on restoring the mass, and tends to zero. The second contribution from the cut is understood 
in the sense of principal value. It is an odd function of /i regular in a vicinity of zero. Consequently, 
it vanishes at /u — )• too. In the fermionic case, the main (quasiclassical) contribution to the 
logarithm of the partition function looks like 

/3J^|1 = ^(1 - 2-'^)(47r)('^+i)/2r {{d + l)/2) Q{d + 1), (113) 

where we have used the basic functional equation for the ^-function. 

To sum over the poles in the case of the Bose-Einstein distribution, it is necessary to isolate the 
contribution from the pole A; = since its contribution diverges in the limit /i — t- 0. This divergence 
is expected as long as we sum over pi in the infinite limits and the contribution to the 0-potential 
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bom Pi = diverges at a vanishing chemical potential. According to formulas (^), (|54|), and (|5 
this pole contributes to the il-potential as 

2c{d) « n''/^^ 



Y,'i^Qr'm2,n'qy\X\) - + j;r(0,|A|(/ - so)), (114) 



d|A|'^/2_ 

where d is the physical space dimension, and we neglect all the terms vanishing at sq — )• 0. The 
contribution from the pole k = is "almost" conformal invariant, i.e., apart from the logarithmic 
term appearing in the expansion of the incomplete gamma function in a neighbourhood of zero, 
its contribution is invariant under a simultaneous dilatation of all the metric components gij. A 
conformal invariant part of the fi-potential does not contribute to the average energy. 

It is useful to assign the parameter A to be equal to —iq^- This choice allows us to terminate both 
the sum over q and the sum over p keeping a few leading terms from these series. A contribution 
from the rest of poles is found in the same way as in the fermionic case considered above. The 
oscillating correction to the fi-potential without the zeroth pole contribution can be cast into the 
form 

s^{{d+l)/2-s) Lis-(rf-i)/2 [e ^) , (115) 
while the quasiclassical contribution reads as 



= ^(47r)('^+i)/2r {{d + l)/2) Q{d + 1), (116) 

where we neglect all the terms disappearing in the limit ^ — t- 0. As a result, in the bosonic case 

= /3f^S + fi^'d + I3nf + c{d). (117) 

All the above mentioned conclusions regarding the oscillating contribution to the fermionic 
potential hold for bosons as well. We also see that the oscillating contribution in the bosonic case is 
more suppressed than in the fermionic case. In both cases, it grows exponentially with decreasing 
sizes of the system and becomes appreciable on the scales of the order of It- The expression for the 
main contribution to the ri-potential both for fermions and bosons coincides with the leading term 



of the high-temperature expansions (pq) and (102). 



Now we take into account that, at zero boundary conditions on the fields, the sums defining the 
partition function are not in the infinite limits, but from unity to infinity. Introducing a notation 
for the sums and recalling a symmetry of the spectrum with respect to the substitution pi — )• —pi 
for any quantum number, we arrive at 

oo oo oo oo 

Pl = l Pd = l A Pai=-00 Pa„=-0O 
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where (ai, . . . ,an), cti = l,d, is a collection of nonequal natural numbers in the increasing order. 
The summation is carried over all such collections including the empty one. There are 2*^ such 
collections in d dimensions. An absence of the sum over some pi in the RHS of (118) means that 
p, = 0. 

Consider separately the zero dimensional contribution entering the RHS of Eq. ( |118D . In the 
case where all pi = 0, we have for bosons 



/30o = -ln(l-e^). 



(119) 



At ^ — 7- 0, the divergent contributions to the bosonic 0-potential are the term (119) and the 
contributions from the second term in Eq. ( |114| ) taken in different dimensions excepting zeroth 
oneQ. It is not difficult to comprehend that a passage to the semi-infinite summation intervals 
reduces for these last contributions to the change of summation limits which also become semi- 
infinite, i.e., from unity to infinity. However, to this end, it is necessary to add to the sum the term 
of such a type taken in zeroth dimension. 



(120) 



because it is not contained in ( 118 ), and subtract it. Collecting this subtrahend term with the 

(-1)' 



zeroth dimension contribution ( 119 ), we obtain 

2-'i(-l)<^+i(r(0, -qlso)/2 + ln(l - e^^)) -^^ ^(te + Inql). 



(121) 



Then, denoting by An := (ai, . . . , a„), by gf^^" the respective block (n x n) of the initial metric and 



so on, we come to the rapidly convergent expansion for the bosonic Jl-potential 



7r("-i)/2 2-'^W2 ^/r(n/2,^2 2 , o 
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+ln„2^ , V ^(Q''^^^' 

Pi---Pd=l A 



'9A, 



.r((n + l)/2)^„ 



2rf+27r(n+3)/2 



CS(<7A„,n+l), (122) 



where the summation over collections A no longer includes the empty one. 

As far as the massless fermions are concerned there are not any problems with divergencies of 
separate terms at zero chemical potential like in the bosonic case. So, we immediately write for the 



The zeroth dimension is excluded since, in this dimension, the whole f2-potential is given by formula (111 
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Figure 1. The ratio of the exact value of the S7-potential taken from its definition to its approximate value given by 



Eq. (L22) for bosons (I) and Eq. (123) for fermions (II). As expected, the asymptotic expansions (122) and (123) 
deviate from the exact values when the condition vRTl) is violated. 



fi-potential 



d—n 



vr 



(n-l)/2 



2^— n 



-(l-2-")r((n + l)/2)C(n + l) 



r((n + l)/2)^„, , V"^'^^'v^ 



2, 



g^2gAj-^r((n + l)/2 + s) 
2'^+i7rs!r((n + l)/2 - s) 



(123) 



where we explicitly single out the contribution from zeroth dimension. Remind that, if the above 
mentioned conditions are fulfilled, all the sums in the obtained expansions can be broken off retaining 
several leading terms. Besides, we can eliminate a summation of identical terms in the sums over 



q in formulas (122) and ( |123| ). It is easy to prove the following combinatorial relation 

oo oo d-n 

2-d^(_l)d-n f^i^qlj = Y, E Y.^''-\-l?Ctnfd-k{qlj, (124) 

A qAn=-°° A qA„=lk=0 

where fn{x) are arbitrary functions. Here we have exploited a metric diagonality. The contribution 
from the empty collection is absent both in the left and RHS of the equality. In particular, if fn{x) 
is independent of n then the RHS reduces to the d- dimensional sum in semi-infinite limits. 

A comparison of the obtained asymptotic expansions for the J7-potential with its exact value is 
presented on Fig. ||. 



2. Nonvanishing chemical potential 



Now we briefly describe how the above results change at the non-zero (electro)chemical potential. 
We shall consider only the fermionic case bearing in mind the electrons in graphene (graphite) . Near 
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the Dirac points, the dispersion law of the electrons and holes looks like |35-| 



e{p) = fl±VF\p\, (125) 

where fl w 0.3 eV, and vp ~ 9.1 x 10^ cm/s. We interpret the first term in the dispersion law 
(|125| ) as the chemical potential. This interpretation is correct only if the form of the dispersion 
law does not change during the studied process (a deformation, in our case). Although we shall 
assume a constancy of the first term in ( |125| ), corrections for the variations of this term can be 
directly included to the final expressions for the O-potential and average charge thereby relaxing 



this assumption. Besides, it was shown in |43, 44| that the dispersion law (|125|) has a gap which 



can be modeled by a small mass entering the dispersion law. This gap also depends on whether the 
graphene lattice is deformed or not. We completely neglect this contribution to the dispersion law 
as long as we are going to consider a graphene specimen with the armchair edges where this gap 
is absent in the tight-binding approximation. Thus, at the temperature /3~^ = 1 K, we have the 
dimensionless chemical potential /i ~ 3.5 x 10^ and the thermal wave length It ~ 2.2 x 10^'^ cm. In 
graphene, there are two branches of the dispersion law (125) and also two spin degrees of freedom 



for each of the branches, i.e., in sum, eight degrees of freedom. For simplicity, we shall study a 
ribbon specimen with sizes (Lx, Ly), a, length of the ribbon Ly being large: Ly ^ It- Then, for the 
armchair boundary conditions p7| , |6l] , 62 1 



Px = ^in + Ak/3), n£Z, A; = {0,1, 2}, (126) 

where k = L^/aQ mod 3 and oq ~ 1.42 x 10^*^ cm is a length of the lattice translation vector. The 
eigenfunctions of the one-particle Hamiltonian implying this quantization condition are represented 
by a composition of two wave functions corresponding to the different Dirac points. Therefore, the 
degeneracy over the spin degrees of freedom only survives. 

The quasiclassical contribution to the fi-potential is expressed through the Hurwitz (^-function 
(see, e.g., |^ ]) and can be cast into the form 



f(S!0 = M(4^)Wl)/2r((<i+l)/2)| J2 

^ V I — n 



.=0 nd + 2-2k) 



-^J' 



d 



Lid+1 {e-^) - 2-^Lirf+i (e-2/^)] |, (127) 



2 

where the spin degrees of freedom are taken into account. The contribution from the cut (7^ is 
independent of the chemical potential and has the form (|110[) multiplied by 



8cos(27rA;/3). (128) 
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The oscillating term is given by 

oo 



(129) 



gl fc=0 

where qm = L^/It and the exponentially suppressed contributions from macroscopic dimensions 
are discarded. Recall that only the poles corresponding to the electron branch of the dispersion law 
(|125| ) contribute to the oscillating part. The holes' poles lie on the "unphysical" sheet. We simplify 
this expression under the assumption that the condition (E^ is fulfilled for the zeroth pole: 



AirqmlfJ' + «vr| » 1. 



(130) 



At the value of the chemical potential presented above and the temperature /3~^ = 1 K, this 
condition implies that we are about to consider the ribbon specimen with a width larger than 
5 X 10~^ cm. Then, making use of the asymptotic expansion of formula (45), we come to 



^ r((rf+l)/2 + n) ^ 



k=0 



[UJk - l^i 



)(rf-l)/2 



{AirqraiuJk - i/i))" 



Li(,+i)/2+n(e-2"'^'"("'=-*^) 



(131) 



where ^ := ^ + k/3qm- The respective contributions to the conserved charge Q are easily derived 
from expressions ( |127[ ) and ( |129| ) with the help of the recurrence relation ([To|). As a result, 



Q2 = ^(4vr)('^+i)/^r((d+ l)/2) {p ^^^^ 



)l-2A;^ 



,d-2k 



\d+l- 2k) 

l + (-l)'^ 



Li,(e-^) -2i-<^Lirf(e-2/^)l , (132) 



and 



Q7 



327r^ ^ r((d-l)/2 + n) ^ 



fc=0 



(133) 



These equations implicitly define the chemical potential jl at the fixed value of the charge Q. Under 
normal conditions (unperturbed crystal), this charge is determined by the above mentioned value 
of the chemical potential /i w 0.3 eV. Of course, there is not any uncompensated charge in the 
crystal under the normal conditions. The non-zero charge Q results from the chemical potential 
redefinition made by us. 

The expansions (|129|) and (|133D are slowly convergent with respect to k when 



27r^g™ < 1. 
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Figure 2. On the left panel: The total chemical potential of electrons in the graphene ribbon of the length 
Ly = 1 cm at = 1 K. The undeformed crystal (normal conditions) is of the width Lx — 20 nm and at 
the chemical potential jl = 0.3 eV. The parameter k entering Eq. ( \12(\ ) is set to fc = (I), fc = 1 (II), and 
k = 2 (III). The small plots depict the corresponding oscillating parts of the chemical potential. On the 
right panel: The total charge Q of electrons in the graphene ribbon of the length Ly ^ 1 cm at the fixed 
chemical potential jl ~ 0.3 eV, the temperature f3^^ = 1 K, and fc = (I), fc = 1 (II), and k = 2 (III). The 
small plots depict the corresponding oscillating parts of the charge. 

For example, at the temperature /3~^ = 1 K and the ribbon width Lx = 10"^ cm, the value of 
this parameter becomes 2Tr^qm w 9.1 x lO"''^. If the condition (pi) takes pace, we may apply the 
Euler-Maclaurin formula to Eq. ( |133| ) to arrive at 



8V9 ^ r((d + l)/2 + n) 

i)/2 2^ — - ^ 



{d+i)/2 n\T((d+l)/2-n 



Li(rf+i)/2+n ( e 



(135) 



(47rgm(7r - ifi))' 

In order to take properly into account the boundary conditions chosen by us, we must divide the 
expressions for the f]-potential and the charge Q on 2^^ = 4. Hence, neglecting the exponentially 



suppressed terms in (132), we finally have 



vr 



^(^2+vrV3), 



2L„ 
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r(3/2 + n) 



■ Im 



(7r-i/i)V2 
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3/2+n I e 



3 — 27r<jm(7r-i/i) 



(136) 



(47rgm(7r - in))' 

There are not the contributions from lower dimensions in the case at hand. Note that the obtained 
formulas hold in the zero-temperature limit /3 — t- 00 and allow simply to find the zero-temperature 
asymptotics of the charge Q. 

Now it is not difficult to describe the basic features of oscillations of the chemical potential. 
Expressing perturbatively from the above equations, we find that the period of oscillations of the 
dimensionless chemical potential with Lx reads as 



'irLxLy 



(137) 



37 



It grows with starting with the value 2lT/fJ'0, where /Uq is the dimensionless chemical potential 
at the normal conditions. This value is twice greater in comparison with the naive estimation. 
Keeping in mind the condition ( |134[ ), we obtain that the amplitude of these oscillations is of the 
order of 




(138) 



It decreases with Lx starting with the value 

C(3/2)(2A^^o)"'^', (139) 

where qm corresponds to the unperturbed crystal. We see from these formulas that the order of 
the amplitude and period of oscillations of the chemical potential /i does not appreciably depend 
on the temperature. The plots of these oscillations are presented on Fig. |2|. The oscillation period 
is rather large at the value of the chemical potential fl given above. So, the oscillations of this 
type cannot be seemingly observed as oscillations and not as the one peak. Inasmuch as the period 
of these oscillations decreases with the width of a ribbon, they are likely to be observed under a 
squeezing of the specimen. 



B. Electrons in a thin metal fllm 



In this section, we consider thermodynamic properties of the conduction electrons in a thin 
metal film. This is a classical problem for the Fermi-liquid theory of metals. It was studied in detail 
both for the quadratic dispersion law (refs |p8|-|30|, |63|| ) and, in the quasiclassical approximation, 
for the dispersion law of an arbitrary form (see, e.g., |^ |9|, |ll-14|). It is well-known that the 



thermodynamic properties of such electrons are described by the Fermi-Dirac distribution with a 
satisfactory accuracy. The dimensionless chemical potential fi entering the distribution should be 
expressed through the number of electrons in the conduction band under the assumption that this 
number is conserved. We shall investigate the metals with the only one conduction band, but all 
the below results are straightforwardly generalized to the case of several conduction bands. Also 
we restrict ourself to the isotropic (in terms of the quasi-momentum) dispersion law and derive 
explicit rather simple expressions for the $7-potential and the number of conduction electrons which 
are valid with an exponential accuracy. This approximation is adequate for the electrons located 
in a sufficiently small vicinity of the conduction band bottom. An analogous but less elaborated 
analysis for a cubic specimen of a metal is given in |S4|. 
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Before we proceed to a calculation of the partition function, we have to make some remarks 
on the way how we are going to take into account a finiteness of the crystal. A presence of a 
surface on the crystal changes appreciably the properties of this crystal near the surface, deforms 
its lattice, and even can lead to a lowering of the two-dimensional symmetry of the bulk specimen. 
A description of these effects from the first principles is quite difficult and, in the most cases, does 
not admit an analytic investigation. Hence, we shall suppose that the defects of the crystal lattice 
appearing due to a presence of the surface do not influence the thermodynamical properties of the 
conduction electrons too much. We shall consider the simplest model of a finite crystal. 

Let V{x.) be the effective potential of the electron for the infinite periodic lattice 

F(x + n%) = V{yi), n' G Z, (140) 

where are the lattice translation vectors. A complete set of solutions to the Schrodinger equation 
for the electron is constituted by the Bloch wave functions 

-|^A + F(x)^ V.p(x) = £.(p)^.p(x), (141) 

where m is the electron mass, s is the band index, p is the quasi-momentum, and £s(p) is the 
dispersion law. The Bloch functions are defined in the standard way 

ipsp (x) = e^P'^/^u^p (x) , tpsp+27rhhi (x) = V'sp (x) , Usp (x + ) = Usp (x) , (142) 

where b* are the basis vectors of the reciprocal lattice: a^iy = . For brevity, we neglect the spin 
properties of the electron and shall take them into account only in counting the degeneracy order 
of energy levels. Now we cut out from the infinite crystal a parallelepiped spanned on the basis 
vectors a^ multiplied by some natural numbers. We shall assume that the effective potential of the 
electron in this finite crystal differs negligibly from F(x). 

In that case, the solutions to the Schrodinger equation obeying the Dirichlet boundary conditions 
can be obtained by a linear combination of the Bloch waves "^^p corresponding to the same energy 
es{p)- It is useful to choose the basis {aj} in the x-space and the basis {b*} in the reciprocal 
p-space. Then this boundary conditions imply the following system of equations on the coefficients 
of the linear combination Csp'. 

^c,pu,p(x) 1^,^0 = 0, J] e^"^^'^/^c,p«,p(x) 1^,^0 = 0, i = T;3, (143) 

spi spi 

where the vectors p lie on the constant-energy surface £s(p) = £s, while rii are natural numbers 
specifying a number of unit cells along each of the basis vectors a^ in the finite crystal. Suppose 
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that the equation, 

es(p) = Es, (144) 

(1 2) 

on the quasi-momentum component in the chosen basis pi has two or less solutions ' (e^, p_l) for 
every i. Together with Eqs. ( |143D , this requirement results in the quantization condition 12| 

\pf = 2Trh—, ki G N. (145) 

In order to satisfy this condition, one need to plot in the p-space a parallelepiped with edges 
directed along the basis vectors b*. The lengths of its edges must be equal to ( |145D multiplied by 
the length of the corresponding basis vector b* and all its vertices should lie on the constant-energy 



surface (144). Of course, it is not possible for any value e<j. If the form of the cut crystal does not 
agree with its periodic structure, the conditions ( |145| ) ought to be understood as a quasiclassical 
approximation to the exact quantization conditions. 

As we have already mentioned, we shall consider the electrons with quadratic and isotropic 
dispersion law in the single conduction band 

.(P) = (146) 

where is an effective electron mass. Also we assume that the reciprocal lattice is cubic. Nev- 
ertheless, all the below formulas are immediately generalized by a redefinition of the constants 
entering them to the case where the constant-energy surfaces of the dispersion law ( |144D have the 
forms of ellipsoid with the principal axes directed along the reciprocal basis vectors. For the given 
dispersion law, the quantization conditions ( |145| ) look in the Cartesian basis like 

P^ = he N, (147) 

where Lj are sizes of the crystal. Passing to the notation used in (|5|), we see that, in our case, the 
metric takes the form 

/ ^^2^2 \ 1/2 

Qij = l^^diagiLl Ll, L^), Id := , (148) 

where Id is the thermal de Broglie wavelength. At the temperature = 1 K and the effective 
electron mass equal to its mass in vacuum, the thermal wavelength is approximately 6.6 x 10~^ 



cm. The roots of Eq. (15) become (the parameter m in Eq. (|T^) is put to zero 



Pfc = {sgn{k)^^'^ + 7r2(2fc + 1)2>^ + ^^2^^2^2k + lY - 1^ , A; G Z, (149) 
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where the sign function in zero is defined as sgn(O) = 1. In accordance with our notation, the 
reciprocal temperature /3 is included to the dimensionless chemical potential At the temperature 
/3^^ = 1 K and the chemical potential jl = 11.7 eV (Aluminium), we have for the dimensionless 
chemical potential /u w 1.36 x 10^. 

First, consider the case where a summation over the quantum number in the logarithm of the 
partition function is carried over the infinite limits (—00,00) . In formula ([74|), we derived the 
quasiclassical contribution for this case 

= -vr^/2^Lirf/2+i(-e'^). (150) 

The oscillating contribution is given by the general formula (p^) . For the nonrelativistic dispersion 
law, the distance between poles \pk+i — Pk\ decreases with an increase of the absolute value of the 
number k, and Imp^ grows. Furthermore, the greater the value of the lesser the distance between 
adjacent poles. That is why we can sum over the poles making use of the Euler-Maclaurin formula 
(P^. In the case at hand, we have two subsequences Sa that may be summed. These are the poles 
Sk with negative and nonpositive k, respectively. Thereby, we come to 

rK-i 

.k=0 

To find the oscillating part of the 0-potential, we need to sum this expression over q. Now we use 
the fact that we consider a thin metal film, viz., we suppose that Ly and are much larger than 
Lx and l^. This assumption considerably simplifies the subsequent formulas for the 0-potential 
and the average number of electrons. In virtue of this assumption, the terms of the sum over q with 
the vectors q^ possessing nonvanishing components along y and z axes are strongly exponentially 
suppressed (see the expansion of the Hankel function (^l])). Then the sum over q reduces eff^ectively 
to a one-dimensional one. 



+ ... (151) 



In principle, we can apply to the sums over q of the terms in Eq. (151) the method developed 



in Sec. |II C 3 , but we shall act in other way. We shall derive a simpler looking expansion under 



the assumption that the condition (43) is fulfilled for q = qm and p = pQ. For example, at the 



temperature /3 ^ = 1 K, the chemical potential fl = 11.7 eV, and the film thickness Lx = W ^ cm, 

47rg™|i3o| ^ 7. (152) 



So, this assumptions is quite reasonable. If q^ is so much small that the condition ( pSj ) is violated 
then the rapidly convergent expansion of the ri-potential can be deduced directly from formula (^ . 
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In that case, we just integrate over macroscopic dimensions and cast out the higher terms of the 
series in the quantum number corresponding to the microscopic dimension. 
Thus, supposing the condition is satisfied and using Eq. (|45|), we get 



^ ^0 V4vr; n!r((d+l)/2-n) 

."-l)/2-n (rf+l)/2-n j , . , g?, 

T- ( '2mqrriPk\ , O-l-i-l-Zfi (^2mgmPK^ 

U+l)/2+n LH^+l)/2+n (e ) + ,^2 J'^+3)/2+n + 1 _ 2n LH^+3)/2+n (e 

A;=o ym 

(153) 

where = Lx/Id- Note that, in odd dimensions, the sum over n terminates and the fulfrhment of 
the condition ( |4^ ) is not necessary. 

Now we should use formula ( |118|) for a passage to the semi-infinite summation limits. As long as 
Ly and Lz are much greater than //j, the lower dimensions give small contributions to the partition 
function. For the quasiclassical contribution, we take into account only the first correction due to 
lower dimensional sums. As for the oscillating term, we retain the leading contribution only. This 



approximation is justified when the condition (43) is satisfied and 



max(l,^)Ly//D > 1, (|po|/gm)^/^V^D > 1. (154) 

The first requirement allows us to discard the contributions from lower dimensions to the quasi- 
classical part of the fi-potential, while the second condition permits to make the same with the 
oscillating terms. 

Then, recovering the spin degrees of freedom of the electrons, we arrive at 

« -^[V^g„Li5/2(-e^) - Li2(-e'^)], (155) 

where S := LyL^- The sum over n in the oscillating contribution terminates for d = 3. So, we 
obtain 

■ K-l 



Pk Li2(e2^^^-Pfc) + ^— Li3(e2^*'?™P*) 



+ Li3(e2-^™) + Lu{e^^'^"^P^) - Li5(e2-^™)|. (156) 

The number K entering this formula should be increased until the value of ceases to change 
considerably. At large dimensionless chemical potentials ^u, as it takes place for the conduction 
electrons in a metal, the Euler-Maclaurin formula is valid with a great accuracy and the parameter 
K can be set to unity or even to zero. The obtained expression is an oscillating function of the 
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size Lx- The periods of oscillations at the fixed dimensionless chemical potential are approximately 
equal to 

ik-h/Repk, k = 0,l,... (157) 

If the imaginary part of the roots pk grows rapidly with k then the contributions with the large 
numbers k are strongly exponentially suppressed. In that case, the oscillations mainly occur at 
the frequencies with k = and k = 1. Of course, this also holds at = 1. At sufficiently low 
temperatures, /i ^ 1, we obtain a physically expected value for the oscillation period 

" (158) 



At the temperature /3^^ = 1 K, the chemical potential jl = 11.7 eV, and = m, the oscillation 
period becomes £ ~ 1.8 x 10^^ cm. 

Now we shall find the number of conduction electrons in a metal. To this end, we need to 
differentiate the obtained expression for the O-potential with respect to ^u. It is easy to do by use 
of the recurrence relation 

= ^4-2^ (159) 
that holds for the nonrelativistic dispersion law. It can be proven by a differentiation of the 



expression (H) with respect to /x and integration by parts (see also for details |64|). Then, under 



the assumptions supposed in obtaining the above expression for the J7-potential, we have 

^° ^ -^[V^'?mLi3/2(-e^) + ln(l + e^)], (160) 
for the main contribution to the average particle number. The oscillating part reads as 

y ln(l - e^^-J^Pfe) - Li2(e2^*«'"P^) ^ Li3(e2'^'^™P^ 



'■D 



(161) 



This formulas allow us to find a dependence of the chemical potential on the sizes of the crystal. 
It is not difficult to express [i = fi{f3,N,L) from Eqs. (|160D and (161) numerically. Though, 



it is impossible to represent the expression for /i(/3, N, L) in terms of known special functions. 
Therefore, we derive here an approximate expression for the chemical potential making certain 
assumptions. First, we suppose that the oscillating part is much lesser than the quasiclassical 
contribution <C ■ Whence 

M°^^-(5^oA^°(M))~'iV°^(/xo), (162) 
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Figure 3. On the left panel: I. The chemical potential of electrons in the thin metal film at the temperature I3~^ — 1 
K and the average particle number A'^ = 1.6 x lO^'' what corresponds to the undeformed metal film with the area 
5=1 cm^, the width Lx ~ 1 nm, and the chemical potential fi — 11.7 eV. The small plots depict the total chemical 
potential at = 10~^ K (II) and jS^-^ = 10^ K (III). The small plot (IV) depicts the quasiclassical part /i° of the 
chemical potential at j3~^ = 1 K. On the right panel: I. The quasiclassical contribution to the average number of 
conduction electrons at /3^^ = 1 K and fj, = 11.7 eV. The total number of conduction electrons in the thin metal 
film with the area 5 = 1 cm^ at the fixed chemical potential /i = 11.7 eV, and the temperatures I3~^ = 1 K (II), 
= 1.5 K (III), = 10"^ K (IV), and /J-^ = 10^ K (V). 



where is the solution to Eq. (160) with respect to fi. In order to obtain an exphcit expression for 
fjP, we also assume that /x ^ 1. In that case, recalling the asymptotic expansion of a polylogarithm 
(see, e.g., 64, ^]), we arrive at 

ttS 



1/2 _ 3 



(163) 



An exact solution to this equation is readily found. However, we give only the approximate value 
of the dimensionless chemical potential 



V 9ttS 



1/3 

+ 6 



(164) 



which is valid when the expression in parenthesis is much greater than unity, i.e., when the film 
thickness is not too much small. Making use of this expression, we have under the same 
assumptions 



r) /V^ 



37rS 
4/2" 



1/3 



(165) 



\'9ttS 

what allows us to obtain the oscillating part of the chemical potential. Now it is a rather simple 
task to find the period of oscillations at a fixed number of the conduction electrons and variable 
chemical potential 



V 8N 



(166) 
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The oscillation period increases with the film thickness. In terms of the oscillation period ( |166D , 
the above used assumption on the film thickness turns into the condition ^ L 

In order to estimate the oscillation amplitude of the chemical potential, we observe that an 
absolute value of the expression ( |161| ) possesses a maximum where the argument of the polyloga- 
rithms entering this expression is close to unity. Despite the fact that the first term in ( |161 ) has a 



logarithmic singularity at this point and the rest two are finite in unity, the main contribution for 
the conduction electrons in metals comes from the second term in Eq. ( [I61D . Under the conditions 
considered by us, this term is approximately 10^ times larger than the first term. The third term 
is much smaller than the second one and gives a vanishing contribution when the argument of the 
polylogarithm is equal to unity. Keeping in mind that the polylogarithm in unity is expressed 
through the (^-function (^), we find the following estimation for the oscillation amplitude of the 
dimensionless chemical potential: 

Thus we see that the oscillation amplitude decreases quadratically with increasing of the film 
thickness. The obtained estimations for the period and amplitude of oscillations are in a good 
agreement with the numerical results presented on Fig. ^. 
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